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"No hay nada más excitante que tratar de entender a la naturaleza, pero por encima de 


todo está la mujer” 


Resumen 


En el año 2004, se descubrió un nuevo material de sólo dos dimensiones con propiedades 
únicas y sorprendentes: el grafeno. Desde su aparición, ha logrado acaparar la atención de 
diversos grupos de investigación, ya que se considera que en un futuro próximo será útil 
en la fabricación de transistores más eficientes que los de silicio, procesadores más veloces, 
nuevos páneles solares, componentes delgados, elásticos y muy resistentes, entre muchos 
otros dispositivos. Visto al microscopio, el grafeno es como una sábana arrugada, es 
decir, plana con ligeras ondulaciones, aunque para grafeno a bajas temperaturas, su 
superficie puede ser considerada como plana. Así que, para el estudio del comportamiento 
de los portadores de carga, se puede emplear la ecuación bidimensional de Dirac. Sin 
embargo, cuando el espacio bidimensional ya no es plano sino ondulado, lo anterior 
deja de ser cierto, pues aquí deben considerarse las partículas en un espacio curvo. Por 
lo tanto, en este trabajo se propone una generalización de la ecuación bidimensional 
de Dirac con el fin de describir los portadores de carga moviéndose en grafeno curvo, 
que es el caso de temperaturas superiores a 10 K debido a la presencia de fonones 
flexurales o grafeno doblado. La interacción se toma en cuenta al considerar una métrica 
inducida con el mismo enfoque que utiliza la relatividad general para la descripción de las 
partículas fermiónicas que se mueven en un espacio-tiempo curvo. La ecuación resultante 
permite incluir de forma natural la presencia de ramas de fonones, así como de un campo 
electromagnético externo. No obstante, para grafeno doblado, modelado mediante un 
perfil geométrico tipo sinusoidal, y que ha sido observado experimentalmente cuando el 
grafeno es crecido sobre un sustrato de hierro, el problema puede ser resuelto a partir 
de una ecuación de Mathieu con un parámetro complejo. Las soluciones resultantes 
muestran el fenómeno de resonancia paramétrica. El patrón de resonancia es la misma 
que en el caso de las soluciones de grafeno bajo radiación electromagnética, que además 
conduce una relación de recurrencia que da lugar a regiones de energías prohibidas, lo 
que nos permitiría controlar y modificar el flujo de los electrones. De esta manera, con 
sólo deformar la superficie de grafeno se podrían construir compuertas electrónicas y 
diseñar nuevos transistores. Los resultados obtenidos se publicaron en la revista Physica 
B: Condensed Matter, bajo el título Bending and flexural phonon scattering: Generalized 


Dirac equation for an electron moving in curved graphene. 
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Capítulo 1 


Introducción 


El grafeno da lugar a propiedades sorprendentes y únicas, tales como alta movilidad 
electrónica [1], alta conductividad térmica [2] y alta transparencia óptica [3], que se cree 
serán importantes para futuras aplicaciones [4-9]. Sin embargo, hay ciertas discrepan- 
cias en los valores de la movilidad electrónica, dependiendo de que si las muestras son 
suspendidas o están en un sustrato [10, 11]. A bajas temperaturas, la dispersión por 
impurezas puede ser la responsable de este efecto, que finalmente conduce a una tran- 
sición metal-aislante, ya que el borde de la movilidad aparece cerca de la energía de 
Fermi [12], como se ha confirmado en el grafeno dopado con H [13]. En contraste, por 
encima 1” > 10K tales discrepancias se cree que son una consecuencia del papel funda- 
mental de las vibraciones de la superficie del grafeno (conocido como fonones flexurales 
[14]), causado por la dispersión de electrones, como se ha demostrado muy reciente- 
mente mediante la aplicación de tensión a las hojas de grafeno [11]. Desde el punto 
de vista microscópico, los resultados de la dispersión generan cambios en las distancias 
entre los átomos, dando lugar a fluctuaciones de los orbitales 7 [15]. Además, dado que 
el grafeno puede ser curvado mediante la aplicación de tensión, se ha propuesto hacer 
ingeniería de tensión con el fin de adecuar y modificar las propiedades electrónicas [16]. 
Cabe mencionar también, que el corrugamiento en el grafeno, junto con la aplicación de 
longitudes de onda larga, produce un pseudo-dispositivo magnético que puede modificar 


la propogación de los electrones [17]. 


Por otro lado, se sabe que el comportamiento de los portadores de carga en el grafeno 
puede ser adecuadamente descrita por la ecuación bidimensional de Dirac (ver apéndice 
A). Dicha ecuación también permite evaluar la interacción entre los portadores de car- 
ga y los fonones, a través de acoplamiento mínimo con un pseudo-potencial [14]. Sin 
embargo, siempre que la descripción de la hoja de grafeno se mantenga estrictamente 


plana, los llamados modos de flexión y de tensión, como resultado de las deformaciones 
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transversales, no pueden ser tomados en cuenta en una forma geométrica natural. Es por 
ello, que en este trabajo de tesis se propone utilizar el conocido formalismo del cálculo 
covariante con el fin de adaptar la ecuación de Dirac a una superficie curva y así aplicarlo 


al grafeno curvo, tomando en consideración la interacción electrón-fonón. 
En los capítulos siguientes se muestra el desarrollo del trabajo de esta tesis. 


El en capítulo 2 se da una breve descripción del átomo de carbono y sus formas alotrópi- 
cas, como lo son: el grafito, el diamante, los fullerenos, los nanotubos de carbano y el 


recientemente descubierto, grafeno. 


En el capítulo 3 se explica qué es el grafeno y que lo hace un material único, así como 
su físca básica (estructura y propiedades electrónicas), método de obtención y futuras 


aplicaciones. 


En el capítulo 4 se estudia la ecuación de Dirac para una superficie curva. Esta nueva 
ecuación ahora se adapta al grafeno curvo con las interacciones electrón-fonón, dando 


como resultado una ecuación de Dirac generalizada. 


En el capítulo 5 se expone una solución particular a la ecuación generalizada de Dirac 
en grafeno curvo, es decir, el grafeno es modelado mediante un perfil geométrico tipo 
sinusoidal, como se ha observado experimentalmente cuando el grafeno es crecido sobre 
un sustrato de hierro. Aquí se llegan a soluciones tipo Mathieu pero con un parámetro 
complejo, el cual nos dice que exite un patrón de resonancia. Este patrón de resonancia 
es característico del grafeno plano sometido a radiación electromagnéticas, que además 


muestra una relación de recurrencia que da lugar a regiones de energías prohibidas. 


Por último en el capítulo 6 se presentan las conclusiones generadas y perspectivas. 


Capítulo 2 


El carbono y sus alotropías 


2.1. El carbono 


El carbono es el elemento químico alrededor del cual ha evolucionado la vida, es el más 
conocido e intrigante de la tabla periódica y la base de toda la química orgánica; todas 
las moléculas biológicas importantes, con excepción de la molécula del agua, contienen 
carbono. Tiene una configuración electrónica 15225?2p?. Esta configuración, de acuerdo 
al modelo de capas, nos dice que hay seis electrones orbitando alrededor del núcleo y 
dos ellos se encuentran en la capa más externa (orbital 2p), ver figura 2.1. Estos elec- 
trones, conocidos como electrones de valencia, presentan la facilidad de formar enlaces, 
lo que le permite al carbono conectarse entre sí de muchas maneras y dar lugar a diver- 
sas estructuras tanto moleculares como cristalinas. Tales estructuras son denominadas 


comúnmente formas alotrópicas del carbono, y aquí se ubican el grafito y el diamante. 





FIGURA 2.1: Configuración electrónica del carbono. 
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2.2. Formas alotrópicas del carbono 


2.2.0.1. Grafito y Diamante 


El grafito, elemento usado para fabricar el lápiz común (Fig. Izquierda 2.2), es la forma 
alotrópica más estable del carbono. Está constituido por láminas planas de átomos de 
carbono ordenadas en forma hexagonal y paralelas entre sí (Fig. Derecha 2.2). Cada una 
de estas láminas están unidas débilmente mediante enlaces de Van der Waals, y esto 
permite a las capas deslizarse horizontalmente con facilidad. Por ello, el grafito posee 
un alto grado de anisotropía, de modo que sus propiedades físicas varían notablemente 
según la dirección en que se haga la medida. Por otro lado, el diamante (Fig. Izquierda 
2.3) es la segunda forma alotrópica más estable del carbono, sin embargo, es un material 
con características físicas sorprendentes, muchas de las cuales derivan del fuerte enlace 
covalente tridimensional entre sus átomos (Fig. Derecha 2.3), que hace al diamante una 


estructura mucho más dura que el grafito. 






Grafito 


Red hexagonal 





FIGURA 2.2: El grafito, elemento usado para lápiz común (izquierda) y láminas planas 
de átomos de carbono que conforman al grafito (derecha). 





FIGURA 2.3: El diamante (izquierda) y su estructura (derecha). 
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Para poder entender cómo es que surgen estas formas alotrópicas del carbono, es nece- 
sario explicar el proceso de hibridación, que consiste en la combinación lineal de orbitales 


atómicos distintos, que dan origen a los orbitales híbridos. 


2.2.0.2. Orbitales híbridos 


De acuerdo con la configuración electrónica del átomo de carbono, se tiene que sólo 
los dos electrones que ocupan el orbital 2p están disponibles para formar enlaces, sin 
embargo, debido a que la diferencia de energía entre el orbital 2s y 2p9 es muy pequeña, 
son posibles distintos tipos de hibridación en el momento en que se produce una in- 
teracción del átomo de carbono con otros átomos cercanos. El proceso de hibridación 
comienza cuando uno de los electrones del orbital 2s salta a un orbital 2p, de tal manera 
que el átomo de carbono se encuentra en un estado excitado (antes de la hibridación). 
Finalmente se produce una mezcla o combinación de orbitales sencillos y orbitales en 
un estado excitado. Dicha mezcla permite formar orbitales híbridos equivalentes (Fig. 
2.4). Dependiendo de la cantidad de orbitales que participen en el proceso de la mezcla, 
existirán varios tipos de hibridación. Estos orbitales híbridos resultantes poseen orienta- 
ciones espaciales bien definidas que dan lugar a fuertes enlaces covalentes llamados 0; es 
decir, los átomos están disponibles para la unión a través de sus electrones en el último 
orbital. 


En la denominada hibridación sp? o tetragonal, que se obtiene al combinar los orbitales 
25, 2P.x, 2Py y 2p2, se forman cuatro orbitales híbridos orientados en el espacio, dispuestos 
en una estructura con forma de tetraedro (Fig. 2.5). Esta hibridación es la base de la 


estructura del diamante. 


Otro tipo de hibridación es la que se conoce como hibridación sp? o trigonal, que se 
obtiene a partir de las combinaciones orbitales 2s, 2p., y 2py, e involucra la formación de 
tres orbitales híbridos en el plano x-y con una separación angular de 120? entre sí (Fig. 
2.6). Esta hibridación da lugar a la estructura hexagonal de los átomos de cada una 
de las capas de las que se forma el grafito. Obsérvese que hay un orbital deslocalizado 
que permance sin hibridar, el 2p,. Este es un orbital perpendicular al plano formado 
por los orbitales híbridos, y se halla disponible para la formación de un enlace 7 con 
otro átomo. Son enlaces covalentes pero más débiles que los enlaces o, ya que existe un 
traslape significativamente menor entre los componentes de los orbitales 2p,, debido a 


la orientacin paralela. 


Por último se tiene la hibridación sp o diagonal, que se obtiene al combinar los orbitales 
2s y 2p, el cual produce 2 orbitales híbridos formando entre sí un ángulo de 180 (Fig. 


2.7). Esto hace que las moléculas que se forman sean lineales. A parte de los orbitales 
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hibridados también se tienen dos orbitales deslocalizados 7 que no entran en el proceso 
de hibridación; ejemplos de hibridación sp los encontramos en la molécula C'Oz. Es 
importante hacer la observación que el proceso de hibridación no se produce en átomos 
aislados, sino más bien, surge de la interacción del átomo de carbono con otros átomos 


cercanos, ya sean del mismo tipo u otros. 





FIGURA 2.4: Hibridación de los orbitales atómicos del átomo de carbono. 


ed 





4 órbitas sp? 


FIGURA 2.5: Orbitales atómicos de valencia 25, 2p,, 29, y 2p2 del carbono (arriba) y 
orbitales atómicos conocidos como híbridos [sp*] = a2s + b2p. + c2py + d2p, donde los 
coeficientes de mezcla a, b, c y d, indican el aporte de cada uno de los orbitales atómicos 
originales a la dirección tetraédrica requerida. Hay 4 posibilidades (abajo). Ref. [18]. 


Capítulo 1 El carbono y sus alotropías 











sp orbital 





FIGURA 2.6: Orbitales atómicos de valencia 2s, 2p,, 2py y 2p2 del carbono (izquierda 
arriba) y orbitales atómicos conocidos como híbridos [sp*] = a2s + b2p. + c2py, donde 
los coeficientes de mezcla a, b y c, indican el aporte de cada uno de los orbitales atómicos 
originales a la dirección trigonal requerida; es decir, se generan tres orbitales híbridos en 
el plano x — y, dirigidos a 120? entre sí (derecha arriba). Obsérvese que aun queda libre 
un orbital atómico de valencia 2p, que no participa en esta hibridación trigonal, por lo 
que efectivamente la forma final es la que se muestra, con el orbital 2p, perpendicular 
a los sp? (abajo). Cada sp? dispone de un electrón para compartir y formar enlaces, ya 
sean con átomos iguales (tipo sigma o) o átomos diferentes, y el 2p, también dispone 
de un electrón para para formar enlaces de características r. Imagen tomada de [18]. 
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FIGURA 2.7: Orbitales atómicos de valencia 2s y 2p.,. del carbono (izquierda arriba) y 
orbitales atómicos conocidos como híbridos [sp] = a2s + b2p.., donde los coeficientes de 
mezcla a y b, indican el aporte de cada uno de los orbitales atómicos originales. Así se 
generan formas lineales de enlaces a 180? entre sí (derecha arriba). Además podemos 
decir que hay 2 orbitales tipo 2p libres que no participan en la hibridación: 2p, y 2p.. 
En la parte de abajo se encuentra la molécula C'Oz, allí se muestra el átomo de carbono 
(hibridación sp más 2 orbitales atómicos originales tipo 2p perpendiculares entre sí) y 
los oxígenos O (hibridación sp? más un orbital atómico original tipo 2p perpendicular), 
uno de ellos girado 90? alrededor de la línea de unión. Imagen tomada de [18]. 


2.2.0.3. Fullerenos y nanotubos de carbono 


A pesar de que el enlace del carbono ha sido uno de los más estudiados, en 1985 se 
descubrió una nueva forma del carbono [19] (de hecho una familia entera de nuevas 
formas). En esta nueva alotropía los átomos de carbono presentan una hibridación in- 
termedia entre el orbital híbrido sp? y sp*. Este tipo de hibridación hace posible que los 
átomos de carbono puedan combinarse formando hexágonos y pentágonos en estructuras 
tridimensionales cerradas. El primer miembro de esta familia y el mejor conocido tiene 
una estructura esférica, compuesta por 60 átomos de carbono (Fig. 2.8). Esta bola de 
dimensión cero y fórmula Céo, similar a un balón de fútbol, se conoce también como 


buckminsterfullereno o simplemente fullereno [19, 20]. 


En 1952 nuevas estructuras del carbono se habían descubierto: “los nanotubos de car- 
bono” (Fig. 2.9). Pero nadie les dio tanta importancia, no fue si no hasta 1991 cuando 
S. lijima [21] los redescubrió para la ciencia occidental y notó una serie de propiedades 
sorprendentes. Los nanotubos de carbono son considerados como una forma alotrópica 
del carbono de dimensión uno. En esta forma alotrópica los átomos de carbono presen- 


tan también una hibridación intermedia. Este tipo de hibridación hace posible que los 
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átomos de carbono puedan combinarse formando hexágonos y heptágonos, dando lugar, 
al contrario de los fullerenos, a una curvatura inversa. Asimismo presenta un orbital 


deslocalizado que permite formar enlaces T. 


Ambas formas alotrópicas del carbono tanto los fullerenos como los nanotubos de car- 
bono han sido extensivamente estudiadas debido a sus interesantes propiedades mecánicas, 
ópticas y eléctricas [22]. Sin embargo, faltaba un compuesto bidimensional del carbono: 


“el grafeno”. 





FIGURA 2.8: Molécula de carbono (C60) conocida también como buckminsterfullereno 
o simplemente fullereno. 





FIGURA 2.9: Nanotubo de carbono. 
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2.2.0.4. Hacia el descubrimiento del grafeno 


Como se ha mencionado anteriormente, ya se conocían varias formas alotrópicas del 
carbono; fullerenos (dimensión cero), nanotubos (dimensión uno), grafito y diamante 
(dimensión tres). Por lo tanto, la búsqueda del grafeno (dimensión dos) era de esper- 
arse, pues el grafito se forma apilando capas planas de átomos de carbono de dimensión 
dos, por lo que cada una de estas capas es precisamente un cristal de grafeno. Afortu- 
nadamente, en el grafito la interacción entre capas es mucho más débil que la interacción 
interna en cada cristal de grafeno, donde los átomos están unidos mediante enlaces cova- 
lentes, el cual es muy fuerte. Esto permite que se puedan utilizar en las puntas de lápiz 
o lubricantes. Entonces podría pensarse que la obtención de grafeno era algo simple de 
realizar. Sin embargo, existían dos problemas: aislar y caracterizar un cristal de grafeno, 
y lograr vencer la supuesta imposibilidad teórica de obtener un material bidimensional 
[23]. Hasta el año 2004 no se creía posible su existencia como entidad aislada, ya que 
se suponía que los cristales estrictamente bidimensionales eran termodinámicamente in- 
estables. El argumento fundamental que reside detrás de tal creencia se debe a L. Landau 
y R. Peierls, ya que setenta años atrás demostraron que tal inestabilidad se debía a las 
fluctuaciones térmicas en redes cristalinas de baja dimensionalidad, lo que produciría 
desplazamientos atómicos comparables a las distancias interatómicas a cualquier temper- 
atura finita [23], dando lugar a la desintegración del cristal. Más tarde Mermin extendería 
este resultado a sistemas bidimensionales más generales, estableciendo que no existe or- 
den de largo alcance en dos dimensiones [23]. Esta hipótesis se refuerza por numerosas 
pruebas experimentales, entre ellas se encuentra el hallazgo de que la temperatura de 
fusión de láminas delgadas decrece rápidamente al disminuir su espesor, provocando 
que la lámina se vuelva inestable para grosores correspondientes aproximadamente a 
una docena de monocapas [23]. Esto fuerza a los cristales bidimensionales a tomar una 
variedad de estructuras tridimensionales que les proporcione estabilidad a costa de la 
pérdida de la bidimensionalidad. De esta manera se suponía que el grafeno únicamente 
podía existir como constituyente básico de otros materiales grafíticos, aunque sin poseer 
una entidad real como objeto aislado. No obstante en 2004, en la Universidad de Manch- 
ester, dos científicos rusos, Andre Geim y Konstantin Novoselov, obtuvieron por primera 
vez láminas individuales de grafeno empleando un método conocido como exfoliación mi- 
cromecánica o “cinta pegante” [24]. Seis años más tarde Geim y Novoselov recibirían el 


premio Nobel de física por su notable descubrimiento. 


Capítulo 3 


Grafeno 


3.1. ¿Qué es el grafeno? 


El grafeno es un material nanométrico de dos dimensiones y el primer ejemplo en el 
mundo real. Consiste en átomos de carbono arreglados en forma de una red hexagonal 
dando la apariencia como la de un panal de abejas (Fig. 3.1). Fue obtenido a partir del 
grafito en 2004 por los científicos rusos Andre Geim y Konstantin Novoselov [24]. Es 
un bloque de construcción básico para los materiales grafíticos. Puede ser envuelto en 
fullerenos (0D), enrollados en nanotubos (1D) o grafito apilado (3D), como se muestra 
en la Fig. 3.2. Posee propiedades eléctricas, térmicas, mecánicas y ópticas jamás antes 
descubiertas por un sólo material, por lo que ofrece y promete una gran cantidad de 


aplicaciones tecnológicas, así como también, mucha investigación en ciencia básica. 





FIGURA 3.1: Estructura hexagonal de los átomos de grafeno (izquierda) y panal de 
abejas (derecha). 
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FIGURA 3.2: El grafeno como bloque de construcción básico. Se puede curvar en 
fullerenos (0D), enrollar en nanotubos (1D) y apilarse en grafito (3D). Ref. [23] 


3.2. ¿Qué hace al grafeno un material único? 


A diferencia de las otras formas del carbono, el grafeno da lugar a propiedades sorpren- 
dentes y únicas, esto se debe a la disposición que tienen sus átomos. Los electrones se 
mueven a través de él como si no tuvieran masa (como los fotones), con una movilidad 
electrónica 100 veces mayor que la del silicio y trescientas veces menor que la velocidad 
de la luz [1, 25]. A temperatura ambiente, posee una enorme longitud de coherencia 
cuántica, es decir los electrones pueden viajar a través de varias micras de grafeno sin 
dispersión; un orden de magnitud mayor que cualquier otro material. A pesar que es el 
material más fino conocido (su grosor no excede el de un átomo de carbón), también es 
el más fuerte jamás registrado: 100 veces más fuerte que el acero [25]. Se ha dicho recien- 
temente que una capa de grafeno podría sostener un elefante sobre un lapicero. Además, 
el grafeno es uno de los mejores conductores del calor (superior que el diamante) [2, 25]. 
En la Fig. 3.3 se muestran las comparaciones gráficas de algunas propiedades del grafeno 


con las de otros materiales. 
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FIGURA 3.3: Grafeno vs otros materiales. Imagen tomada de [25]. 


Aparte de tener alta movilidad electrónica [1, 25], alta conductivad térmica [2, 25] y ser 
muy resistente [25], el grafeno es aún más atractivo debido a su transparencia óptica: 
“sólo el 2.3% de la luz blanca es absorbida” [3]. Es tan denso, que ni siquiera el gas helio 
(el átomo más pequeño) lo puede atravesar [26]. Se extiende con facilidad y es muy flex- 
ible, lo que permite la capacidad de doblarse (otros materiales resistentes se romperían). 
Asimismo es muy sensible a cualquier molécula que se deposite en su superficie [27]. 
También si se le aplica una señal eléctrica de cierta frecuencia, genera otra onda del 
doble o el triple de frecuencia, es decir, funciona como un multiplicador de frecuencias; 


por lo que permitirá trabajar a frecuencias más altas de las actuales, en el rango de 500 


a 1000 gigahertz [28]. 
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3.3. Física básica del grafeno 


3.3.1. Estructura 


La estructura cristalina del grafeno consiste de una red hexagonal. Esta red puede verse 
como la composición de dos subredes unitarias triangulares no equivalentes A y B, con 
una base de dos átomos de carbono por cada subred, como se muestra en la Fig. 3.4 
(izquierda). La utilización de estas dos subredes triangulares proviene del hecho que la 
red triangular es una red de Bravais, mientras que la red hexagonal original del grafeno 


no lo es. 


Los vectores unitarios de cada red triangular se pueden escribir como: 
a a 
a=3(8,V3), a =>3(3,-v3), (3.1) 


donde ay = 1.42 4, es la distancia carbono-carbono. 


pr 
/ xA 
e 





FIGURA 3.4: (Izquierda) Los átomos azules forman la subred unitaria triangular A y 
los átomos amarillos la subred unitaria triangular B. (Derecha) Los puntos de Dirac K 
y K' de la zona de Brillouin del grafeno. Imagen tomada de [14]. 


Por otro lado, los vectores unitarios de la red recíproca son: 


21 


br a. des =—(1,V3) (3.2) 


340 a0 


Los tres vectores de los átomos vecinos en el espacio real están dados por: 


di 3 (1, V3), d= (1, -V3) y 03 =-—ao(1,0) (3.3) 
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Existen también los puntos K y K”, como se muestra en la Fig. 3.4 (derecha); llamados 


puntos de Dirac, que se encuentran en las esquinas de la zona de Brillouin del grafeno: 


27 27 27 27 


a . Hess l 3.4 
30 3v3a0 Ga 3vV3a0 ) 








K =( 


Todas las propiedades eléctricas que dan al grafeno ese carácter especial se centran princi- 
palmente en torno a estos puntos, ya que su particular estructura de bandas electrónicas: 
el famoso cono de Dirac, formado por la banda de valencia y la banda de conducción, 
se juntan en estos puntos a una energía conocida como energía de Fermi. Obsérvese [29] 
(Fig. 3.5) que si el nivel de Fermi (ef) está en la banda de conducción nos encontramos 
ante un metal, ya que los electrones circulan libremente por esta banda. En el caso de los 
semiconductores y los aislantes el nivel de Fermi está entre ambas bandas de conducción, 
y únicamente se diferencian en la anchura de la brecha, llamada gap o banda prohibida. 
En los semiconductores tenemos un gap no demasiado grande, lo que permite que si se 
les da la suficiente energía puedan conducir la corriente eléctrica; mientras que en los 
aislantes, donde el gap es mucho mayor, es necesario un aporte energético superior que 


dificulta la conducción. 


GRAFENO SEMICONDUCTOR AISLANTE 


METAL 


ENERGIA 





FIGURA 3.5: Bandas de energías. Metal, grafeno, semiconductores y aislante. Imagen 
tomada de [29]. 


3.3.2. Propiedades electrónicas 


El grafeno, al igual que el grafito consiste en la hibridación sp?. Dado que tres de los 
electrones de carbono se usan para el amarre covalente gd en el plano de la red hexagonal 
(Fig. 3.6), el cual da la solidez a la estructura reticular, se tiene que el responsable de 
las propiedades electrónicas es el electrón débilmente ligado en el enlace 7. (Fig. 3.7). 


Mediante un hamiltoniano de enlace fuerte (tight binding) aplicado al electrón T, se 
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FIGURA 3.6: Los enlaces fuertes o de la hibridación sp? determinan la estructura de la 
red hexagonal del grafeno (tipo panal de abeja) y le da solidez. 





FIGURA 3.7: El enlace 7 determina las propiedades de transporte del grafeno. 


llega a que la relación entre la energía (E) y vector momento K es [14]: 


Ex(K) = +ty3+J(K) — 1 f(K), 
con f(K)=2cos(K - ar) + 2cos(K - az) + 2cos(K - (a2 — a7)). 





(3.5) 


El signo “+” en esta ecuación se refiere a la banda superior (o banda de conducción), 
conocida como la banda T y el signo “—” indica la banda inferior (o banda de valencia), 
conocida como la banda 7*; de la estructura de banda del grafeno. El parámetro t = 
7 =3eV [30] es la energía de salto al primer vecino y el parámetro 1 es la energía de 
salto al segundo vecino. La figura 3.8 muestra la estructura de banda completa para el 
grafeno; en esta misma figura se puede ver que las bandas se tocan en los puntos de 
Dirac, exactamente a la energía de Fermi, que en este caso coincide con la energía cero. 
Notemos también que cuando t' es cero, el espectro es simétrico alrededor de la energía 
de Fermi, por esta razón se considera al grafeno como un semiconductor con banda 
prohibida de energía de tamaño cero, mientras que para un valor t' finito la situación 


cambia debido a la ruptura de esta simetría. 
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FIGURA 3.8: estructura de banda completa para el grafeno. Imagen tomada de [14]. 


Desarrollando en serie la Ec. 3.5, cerca de los puntos de Dirac como q = K + k, con 


| k [<c] K |, se tiene que [14]: 


E¿(k) =+0p | k | +e((1k]/|K |), (3.6) 





donde k es el momento medido relativo a los puntos de Dirac y vp representa la velocidad 
de Fermi, que está dada por vp = 3ta/2 = vp = 1x100 m/s. La velocidad de Fermi en 
estos puntos de Dirac es una constante que no depende de la energía o el momento, como 
es típico de los semiconductores con curvas parabólicas en su relación de dispersión. Este 


resultado fue obtenido por primera vez por Wallace en 1946. 


Como las propiedades electrónicas dependen de los electrones con energía cercana a 
la energía de Fermi, se puede despreciar t' y basta solamente estudiar sus entornos. 
Allí el grafeno depara una sorpresa cuando k es muy pequeño. Por lo que, únicamente 
la aproximación a primer orden en k es importante. De esta manera, considerando tal 
aproximación, obtenemos una dispersión lineal en la energía, el cual es característica 
de la electrodinámica cuántica QED para fermiones de Dirac sin masa, excepto por el 
hecho de que, en el grafeno, los electrones ahora tienen una masa efectiva cero y se 
mueven con una velocidad vp, que es 300 veces menos que la velocidad de la luz (c) 
en el vacío [1, 25]. Por lo tanto, muchas de las propiedades de la QED pueden aparecer 
en el grafeno, pero a velocidades más pequeñas. Una de estas propiedades que aparecen 
es el efecto Hall cuántico anómalo [23], el cual ha sido comprobado experimentalmente 
por diversos grupos [1, 31]. Esto es tal vez, la demostración más interesante de que los 
electrones en el grafeno se comportan como fermiones de Dirac sin masa. Además, los 
electrones pueden ser transmitidos con probabilidad 1 a través de una región prohibida 


[32], lo que se conoce como paradoja de Klein y que aparace también en la QED. 
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3.4. Método de obtención 


El grafeno en estado libre fue obtenido por vez primera en la Universidad de Manchester, 
en 2004 por los científicos rusos A. K. Geim y K. Novoselov [24], mediante la técnica 
conocida exfoliación micromecánica o cinta pegante, consiste en la separación de la 
capa más externa de un sólido en laminillas. Por lo tanto, la exfoliación micromecánica 
utilizada para separar capas de grafeno se basa en tomar una muestra de grafito con 
cinta adhesiva; cabe mencionar que un milímetro de grafito consta en realidad de tres 
millones de capas de grafeno. Esta cinta adhesiva con muestra se pega y despega. En 
principio se obtienen muchas capas unidas de grafeno, pero cuando se repite el método 
de pegar y despegar unas 20 veces, se encuentran capas aisladas de grafeno del grosor 
de un sólo átomo de carbono (Fig. 3.9). Para corraborarlo, las muestras se fijan a una 
placa de silicio oxidado y se analiza mediante microscopía electrónica de tunelamiento, 


microscopía electrónica de barrido, microscopía de alta resolución y efecto túnel [33]. 





FIGURA 3.9: Exfoliación micromecánica, método para obtener grafeno. 
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3.5. Futuras aplicaciones 


Las propiedades que exhibe el grafeno ha despertado el enorme interés de investigadores 
y empresarios por producirlo a escala industrial dadas sus potencias promisorias en 
aplicaciones tecnológicas. La investigación en este material se ha centrado en las futuras 
aplicaciones en nano-dispositivos [4-6]. Algunos artículos publicados han sido acerca 
de transistores, fotodetectores y sensores [7-9]. De todas las aplicaciones sugeridas de 
grafeno, el uso de grafeno como un electrodo delgado parace el más cercano a emerger. 
Ello se debe a las propiedades eléctricas (es un excelente conductor de la electricidad) y a 
las propiedades magníficas en la región de la luz visible, con una transparencia mayor al 
90% [3]. Los electrodos transparentes se requieren en una gran variedad de aplicaciones, 


tales como la pantalla táctil y pantallas de cristal líquido. 


Capítulo 4 


Ecuación de Dirac generalizada 


en grafeno curvo 


Dado que los portadores de carga en el grafeno plano son descritos por fermiones de Dirac 
sin masa (ver apéndice A), es natural preguntarse si es posible modificar la ecuación de 
Dirac, teniendo en cuenta la interacción entre los fonones flexurales. Dos caminos se 
pueden seguir con el fin de responder a esta pregunta. Uno de ellos es el enfoque usual 
del enlace fuerte (tight binding) que utiliza el desarrollo de Taylor de la integral de 
traslape en el campo de desplazamientos [14, 15, 34]. Aquí se presenta un punto de vista 
alternativo, en el que se incluye la interacción entre los fonones flexurales, haciendo la 
observación de que una hoja de grafeno puede ser considerado, como un espacio curvo. 
La ecuación efectiva debe ser covariante; es decir, invariante bajo transformaciones de 
Lorentz. La ecuación deseada se puede considerar como la ecuación de Dirac en espacio- 
tiempo curvo, como en relatividad general. Además, dado que los modos de vibración 
son descritos por un potencial vectorial [14, 15, 34], nuestro enfoque permite describir 
todas las ramas de fonones y el potencial electromagnético, en una misma ecuación. Hay 
que tener en cuenta que el problema de la ecuación bidimensional de Dirac incluyendo un 
potencial vectorial electromagnético se ha resuelto recientemente [35, 36]. Como veremos 
al término de este capítulo, el enfoque actual tiene ciertas ventajas sobre el desarrollo 


de Taylor de los parámetros de enlace fuerte. 


A continuación vamos a esbozar brevemente las ideas detrás del enfoque actual. Para 
erafeno a bajas temperaturas, su superficie puede ser considerada como plana (Fig. Izq. 
4.1), por lo que la ecuación de Dirac no modificada es utilizada para describir la evolución 
de las cargas en la hoja de grafeno. Esta se escribe como sigue: 


19) de % 
hh Y =HVY con H =-—ihuplyV y + y2V yl, (4.1) 


20 
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donde vf es la velocidad de Fermi y 





FIGURA 4.1: Grafeno plano (temperaturas bajas) y grafeno curvo (T” > 10K). 


Las matrices anteriores muestran en este momento que no hay ninguna diferencia entre 


los índices covariantes y contravariantes, debido a que: 
yn + ye = 0305 + 050; = 2951, (i,j =x,4), (4.2) 


con gi; = 9;;, por lo que, obviamente, la métrica contravariante elevando los índices es 
también gY = 5%, Pero esto deja de ser cierto cuando el espacio bidimensional no es 
plano, sino ondulado (Fig. Der. 4.1), como es el caso para grafeno a T > 10K'. Esta es 
la razón por el cual la ecuación de Dirac en dos dimensiones debe ser ahora generalizada 
con el fin de incorporar el hecho de que la métrica en la superficie de restricción no es 
más plana, sino curva. Afortunadamente, el problema de la formulación covariante de 
la ecuación de Dirac en un espacio curvo ha sido profundamente investigado desde hace 
mucho tiempo [37, 38], de manera que podemos seguir los mismos pasos en este caso 
concreto: introducir la métrica no euclidiana de dos dimensiones en la hoja de grafeno, 
y luego adaptar el álgebra de Clifford y encontrar la conexión de Christoffel, y montar 


todo esto en la versión covariante de la ecuación de Dirac: 


HU A [FV + 9V]V. (4.3) 


Aquí no sólo la métrica contravariante se deforma, sino también las matrices-y deben 
ser modificadas a fin de satisfacer las nuevas relaciones anticonmutativas con la métrica 
inducida, en lugar de la métrica plana. Además, V;¿ no sólo contiene la interacción 


electromagnética y el plano de fonones visualizado por el potencial vectorial, incluido de 
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manera usual en la invarianza de norma, sino también los símbolos de Christoffel de la 


métrica Ji: 


V¡v = (0; A v 


= AU+D7g P v, (4.4) 


donde 7; = 0; — eA; y los símbolos de Christoffel se definen como de costumbre, por 


medio de la métrica modificada: 
Ei das > > 
je = 99 (Ojdmk + Ox Gjm — OmIjr) (4.5) 


y .. es el valor matricial del tensor antisimétrico definido por medio de las matrices- 


gamas modificadas: 


Y. = ¿Bn (1.6) 


Este término a menudo se llama ”conexión espinorial” [37, 38]. 


Ahora, como se demuestra en el apéndice B, usando los términos hasta segundo orden 


en las potencias de la deformación 0;f, el resultado es como sigue: 


ñ = UF y fa — Uy He (O, f) Te | AO 0 
ii j B (Oy F)Tz — (Oy f )ñy 


—0p(ÚS-GNVS 4) — ihr [Eo VOLS) +(0,8 7] - 418 - IAS, (4.7) 
donde Y = (9,,0,) y Ti = 05 — eAj (j = x,y) es el momento generalizado. 


Aquí podemos empezar a investigar las soluciones que describen el comportamiento de 
los portadores de carga en la hoja de grafeno deforme, interacturando con los fonones 


flexurales. 


Para mostrar cómo se puede utilizar la Ec. 4.7, consideremos que el campo de desplaza- 
mientos puede ser escrito usando un conjunto simple de funciones base proporcionada 


por las ondas estacionarias [39]: 


1-5 DE eu) 10 0%, cos( (7 - E=oÑa 0+Q sind > H= w (q )t)], (4.8) 
1,9>0 


donde 87) es el vector de polarización para un vector de onda 7, y es la rama 
fonónica, wulg) la relación de dispersión y R es la posición. La notación q > 0 indica 


que qz > 0 y qy > 0. 07 y aq son los operadores dados en términos de los operadores 
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fonónicos de creación (aniquilación) [39] a,,, ¿(al dd 7) 
(0) _ 1 (2) ¿0 
44 V2Mw (4) (al +a( Y) (1) 


donde q: corre sobre c y s, y M es la masa del átomo de carbono. Por lo tanto, para 


fonones flexurales la función f(x, y,t) está dada por: 


fe at= Ss $ [06 ¿cosóy + Q (2, sing] y (4.10) 

q>0 
donde y = F significa que estamos trantando con fonones flexurales, y la fase (y se 
define como $, = 7 * R dq): Para la rama flexural, wp(4) = ap || 7 [1? en donde 
ap = 4,6x107m?/s. ao, y 00 tienen unidades de longitud, por lo que f(x, y, t) 
tiene las mismas unidades, mientras que 09, f y 0, f son adimensionales. A primer orden 


la ecuación de Dirac generalizada sólo necesita como entrada las derivadas parciales de 


Fiz, Y, t): 


Ox f = G y qx [Q ¿cosóa + Ql 6), sindy] (4.11) 
q>0 

di G y dy [Q ¿cosóa + Ql 6), sindy] (4.12) 
q>0 


Por ejemplo, si tenemos sólo la corrección lineal de la Ec. 4.7 y sustituimos las Ecs. 4.11 


y 4.12 , el hamiltoniano es el siguiente: 
S 0 Ta — My y E 
H=vp| . Ml 2 [06 cosóy + Q (2 singy|x 


Ta + Uy 0 
QaTo + qyTy 0 ) (4 13) 


0 — Qu Ta NN qyTy 


Usando el hecho de que la velocidad de Fermi es mucho mayor que las velocidades de los 
modos flexurales, podemos despreciar la dependencia del tiempo y aplicar la teoría de 
perturbaciones para resolver la ecuación. Para entender la naturaleza de las soluciones, 
en la siguiente sección vamos a considerar el caso particular de una superficie con perfil 


geométrico tipo sinusoidal. 


Capítulo 5 


Soluciones aproximadas 


Como se dijo en el capítulo anterior, aquí estudiaremos una de las superficies curvas más 


simple, dado por una forma geométrica sinusoidal: 


f(x, y) = acos(q - r), (5.1) 


donde a es la amplitud de la deformación. Este corrugamiento ha sido experimentalmente 


observado en grafeno crecido sobre un sustrato de hierro (110) [40]. 


Al introducir la función deformación f en la ecuación de Dirac generalizada a primer 


orden (Ec. 4.13), se obtiene: 


(T, — iy) YB + acos D¿(q4Ta + qyTy) Va = (ih/vp)O DA, 
(To + Ty) VA + QcOs Py¿(—QrTe == QyTy) VB = (ih/vr)0VB. (5.2) 


En el caso particular en que no hay dependencia del tiempo en f(x, y;t), la evolución 
espacial y temporal puede ser separada. La parte temporal de la función de onda tiene 
la forma exp(—iwt), donde w es una frecuencia. Por lo tanto, sólo nos quedamos con la 


parte espacial de la función deformación a cos(q : r). Eligiendo q = (q, 0), se obtiene: 


(Ox — 10y)VB + acos(q1)90VA = EVA, 
(0, + 59,)w4 — acos(qu)a0.18 = Epa, (53) 


donde el parámetro E se define como E = ihw/hvp = ¿e/hvp, donde e es la energía, 
y ahora YA y Yg denotan sólo la parte espacial de la función de onda. Con el fin de 


separar las variables, supongamos que la solución estacionaria es de la forma 
Valz,y) =F(y)Da(o), va = Fly) plz), (5.4) 
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donde F es una función regular en el origen. Entonces las Ecs. 5.3 toman la forma: 


F(y)g(x) -¿F(y)0g(x) + aqF (y) cos(q1)0', = EF(y)Da, 
F(y)0'4(2) +5F (y) 4(1) — aqF (y) cos(q1) Dg = EF(y)0 p. (5.5) 


Después de dividir por F(y), vemos que las variables se separan fácilmente si ponemos 


E = Const, 
lo que significa que F(y) es una función exponencial con exponente real o imaginario 
puro, dependiendo del signo de la constante. Elegimos el signo negativo con el fin de 


recuperar el caso de grafeno plano en el límite apropiado, lo que lleva al siguiente ansatz: 

















Valor, y) =D (o), pala, y) = Ola). (5.6) 





La función 4 tiene dos signos, ya que, en principio, Va y Yg son las componentes 
de un espinor, y tenemos que construir dos soluciones, una con la energía E y la otra 
con —E. Para el caso de grafeno plano, se sabe que uno puede pasar de una solución a 
otra, cambiando el signo de una componente, debido a la simetría de la red. Estas dos 
soluciones son necesarias para representar electrones y huecos. Aquí nos concentraremos 
en electrones, aunque la solución para los huecos se puede encontrar en una forma similar. 


Por esa razón, en lo que sigue vamos a buscar la solución para vi. La inserción del ansatz 


5.6 en la Ec. 5.3 da como resultado: 


(0, — (¿K,y)Pg(x) + qa cos(qu)0,Pa(x) = Eba(x), 


(O, + i(1K,y)b4(x) — qacos(q1)0,Pg(x) = Ebg(zx), (5.7) 

lo que equivale a 
[O, + Ky]Pz(x) + qa cos(q1)0,Pa(3) = Ebalz), (5.8) 
(0, — Ky]Bg(x) + qacos(q1)0,L (3) = EbaA(z). (5.9) 


Este sistema de ecuaciones es lineal, por lo que podemos obtener la solución general 
sumando soluciones parciales con Ky dado, es decir, utilizando el desarrollo de Fourier. 
Para tener una idea elemental de las propiedades de las soluciones, vamos a examinar 
en primer lugar el modo cero, que corresponde a una corriente incidente en la dirección 
x, dejando Ky = 0. Para otros modos de transporte (con Ky X 0) la diagonalización 
no es sencilla y sólo se puede lograr por aproximaciones sucesivas. En este trabajo sólo 
tomamos en cuenta el caso de Ky = 0, que es la más pertinente en relación con un 


experimento de conductividad típica. Por lo tanto, con el modo cero, la diagonalización 
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del sistema 5.9 es especialmente fácil, ya que toma la forma simplificada: 


09PB(1) + qacos(q1)0,Pr(x) = ED (2), 
0D 1(1) — qa cos(q1)0,Pp(1) = Edp(z). (5.10) 


Este sistema se diagonaliza por derivación de una de sus partes con respecto a x y susti- 
tución posterior en la segunda ecuación. El resultado ordinario de la ecuación diferencial 
de segundo orden es común para ambas componentes Pz y Pp: 


2 do 
1 + a2q? cos*(qu 0 2a* q? sin(qx) cos(qu == aqsin(qu)] BA (5.11) 
dx? dx 





o, en una forma más estándar, 


db, alg? sin(2q1) db,  E[E—aq*sin(qu)) 
dar [1 + a2q? cos?(q1)] de [1 + a2q? cos? (qx) 





dA=0 (5.12) 





y la misma ecuación para 9 pg, pero con el signo + en el tercer término. A continuación 
introducimos un nuevo parámetro y = aq?, que da lugar a una energía efectiva, el cual 


puede entenderse a partir de la relación 
[e=£=zZtan, (5.13) 


donde e es la energía y p es el momento. Este parámetro resultar ser de utilidad teniendo 
en cuenta diversos límites físicos de la energía, como se verá más adelante. Podemos ahora 
reescribir las dos ecuaciones como: 


dba aqnsin(2qx) dd,  E[E-—nsin(qu)] 
dx? [L+ancos?(q1) de  [1+ancos?(qx)) 





dA=0 (5.14) 


dd y aqnsin(2q1) dbg  E|E=+mnsin(qz)) 
dx? [1 + an cos?(qx)] de [1 + an cos?(qx1)) 





dg=0 (5.15) 


Ahora recordemos que en nuestra notación, el parámetro E se ha convertido en un 


imaginario puro: 
E=iT—=i—, (5.16) 


donde e es la desviación de energía del cono de Dirac (ep), es decir, la energía total de 
electrones €s €totar = € + €p. Estas observaciones son importantes a tener en cuenta para 


comprender la naturaleza de las soluciones en lo que sigue. 
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Notemos que la Ec. 5.12 es de la forma estándar 
y" + P(u)y' + Q(2)y =0. (5.17) 
La solución 
y(=) = Z(a)e12S Plodaz (5.18) 
reduce la ecuación a una forma más simple: 
Z"+C(2)Z =0, (5.19) 
donde 
Cía) =Qla) — ¿P (2) - [P2(a) (5.20) 


Con el fin de tener profundización respecto a las propiedades de las soluciones, vamos a 


tratar el caso de una deformación pequeña en el grafeno, es decir, una amplitud a pequeña 


de la superficie f. Tomando esta consideración se puede prescindir de las potencias de 


a superiores o iguales a dos (an < 1). 


Suponiendo que el parámetro E es del mismo orden que el parámetro an (ya que cerca 


del cono de Dirac e << 1), se obtiene la versión simplificada linealizada de nuestras 


ecuaciones: 


29 do 
E aqn sin(2qu) 2 — E(E — nsin(qr))B4 =0. 
de? dx 





Ahora tenemos 


P(=) =—qansin(2q1),  Q(u) =-E(E— ysin(qu)) 


fro = =an | sin(2anjda = a cos(2q1) 


de modo que 


Da(u) — Z(a)e n/4cos(Qaz) 


La función Z(x) satisface la siguiente ecuación diferencial lineal ordinaria: 


3 


Z" + | Ensin(qu) — E? — y? cos(2qx) — 53 sin?(2q1)| Z =0. 


(5.21) 


(5.22) 


(5.23) 


(5.24) 


(5.25) 
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y la inserción de la expresión 5.16, se obtiene 





" € 3 l€ a 2 an? , 
LA nsin(qx) — ny“ cos(2q1) —- — sin(2q1)| Z =0. (5.26) 
hup hup 4 





Fácilmente se comprueba que en el límite de la geometría plana, cuando no hay defor- 


mación (a = 0), la ecuación para Z es: 


2 

€ 

Z" + (53) 2=0, (5.27) 
huvp 

La solución de la encuación anterior representa una onda monocromática que se propaga 


en la dirección x, 
$, = e wt-(e/hvp)x) e ¿ui Kez). (5.28) 


es decir, es la solución del problema de las partculas libres (como era de suponerse) y 


donde se ha puesto 
€ 


Ka=— 
z hu? 


(5.29) 


el cual tiene dimensión, cm”!. 


Despreciando ahora el término cúbico 1?, podemos escribir la siguiente ecuación aprox- 


imada para Z(x): 
2" +|[K? +1iKoysin(qu) — 1” cos(2q1)) Z=0 (5.30) 


Dos límites que son de especial interés pueden ser tratados por separado y dar lugar a 


conocidas ecuaciones: 


I) Cuando K, < y, correspondiente a los electrones entre la energía de Dirac y la 
energía eficaz an (vale la pena recordar que cerca del punto de Dirac, K es casi cero, ya 
que mide la desviación del impulso Kp en cada cono de Dirac [14], con lo que el impulso 


completo es hKrotar = R(Kp + K). Entonces la Ec. 5.30 se reduce a 
Z" — y? cos(2q1)Z =0 (5.31) 


que es inestable por naturaleza, debido a que no hay término constante respecto de la 
frecuencia fundamental en la ecuación de Mathieu: aquí el término oscilante cambia de 
signo y la solución de Z no se puede mantener acotada. Esto significa que los electrones 
cuya energía es demasiado baja (K, < y) no son capaces de propagarse en la hoja de 
grafeno. Notemos también que la inclusión del término lineal y conduce básicamente al 


mismo tipo de solución no acotada. 
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II) Cuando K., > 1, que corresponde a energías altas (más de la energía de Fermi, 


Ky > Kp), la ecuación para Z se convierte en 
Lie Pl + ia siga) Z=i (5.32) 
T 


Esta ecuación tiene la forma de una ecuación de Mathieu con un parámetro complejo, 
que aparece por ejemplo en el problema de corrientes de Foucault en un solenoide elíptico 


[41]. Definiendo una nueva variable de tal manera que 
qu = 2% + (7/2). (5.33) 
La Ec. 5.32 puede reescribirse como 
Z" + [A +28 cos(2q3)] Z=0 (5.34) 


donde los parámetros A y s son tal que: 








4K? 2K 
A= a $= sn 
q q 


(5.35) 


es decir, hay una condición suplementaria en la ecuación de Mathieu, ya que A(s) = 
(q/1m)?s?. Las soluciones de la Ec. 5.34 están dadas por las funciones complejas de Math- 
ieu [41] cer(z, —is), cei(t, —is), ser(E, —is), y sei(í, —¿5). Por lo tanto, la solución puede 


escribirse como una superposición de soluciones con coeficientes C¿ y Dj, 
Z(2) = > ¿[C;ces(%, —is) + Djsej(£, —15)), 
donde, por ejemplo [41], 
00 
cesnlE, —is) = ceraní[T, —is) + iceionl[ E, —is) = 9 Aem) cos(2r2,—is) (5.36) 
r=0 


2n hi Ñ 2 e k A a 
y AR ) es un coeficiente complejo, con una relación de recurrencia definida a partir de 


una variable auxiliar Vo,., 


Año 
Ae 65) (5.37) 





de tal manera que 


(dr? — Ap)Var—2 + s(VarVar-2 +1) =0. (5.38) 
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Las soluciones cen +1(1) tienen una naturaleza similar. Un segundo conjunto de solu- 


ciones son obtenidos de la forma [41] 


00 
seonlT, —is) = seranl[i, —is) + iseignl[T, —1s) = Y Ben sin(2r2,—is) (5.39) 
r=0 


Las soluciones resultantes muestran el fenómeno de resonancias paramétricas. El patrón 
de resonancia es la misma que en el caso de grafeno bajo radiación electromagnética [36]. 
Cabe mencionar que la resonancia paramétrica es un fenómeno que tiene que ver con 
las vibraciones debido a excitaciones externas, en donde uno o más parámetros que lo 
definen son variables en el tiempo. Cuando la excitación es aproximadamente dos veces la 
frecuencia natural de vibración, la oscilaciones tienden a incrementar y por consiguiente a 
producir inestabilidad en el sistema. Entonces, para el caso del grafeno, las excitaciones 
de los átomos de carbono, debido al incremento de la temperatura o por medio de 
las deformaciones transversales, podrían ser los responsables de grandes respuestas en 
el sistema. Esto es de suma importancia en el área de ingeniería mecánica, ya que 
conociendo los requerimientos de un dispositivo, se pueden establecer consideraciones de 
diseño en base a las excitaciones paramétricas. Por otro lado, por medio de la relación 
de recurrencia 5.38, pueden aparecer regiones de energías prohibidas [36] que podría 
permitirnos controlar y modificar el flujo de los electrones, y ser útil para construir 
compuertas electrónicas. Esto es de gran importancia para la fabricación de nuevos 


transitores a base de grafeno. 


Capítulo 6 


Conclusión 


Se dio una breve descripción del átomo de carbono y como a partir de este elemento 
surgen las diferentes estructuras cristalinas: el grafito, el diamante, los fullerenos, los 
nanotubos de carbono y el grafeno. Esta última estructura es fuente de motivación de 
muchos artículos de investigación y por si fuera poco del premio Nobel de Física 2010; 
debido a sus sorprendentes propiedades eléctricas, térmicas, mecánicas y ópticas. Tiene 
además una característica peculiar de la QED, los portadores de carga en el grafeno a 
baja energía, obedecen la ecuación de Dirac; es decir, los electrones se mueven a través 
de él como si no tuvieran masa (como los fotones), con una movilidad electrónica 100 
veces mayor que la del silicio y trescientas veces menor que la velocidad de la luz [1, 25]. 
Sin embargo, hay ciertas discrepancias en los valores de la movilidad electrónica, de- 
pendiendo de que si las muestras de grafeno son suspendidas o están en un sustrato 
[10, 11]. Por encima T' > 10K se cree que tales discrepancias son una consecuencia del 
papel fundamental de las vibraciones de la superficie del grafeno (conocido como fonones 
flexurales [14]); causado por la dispersión de electrones. Por otro lado, el grafeno puede 
ser curvado mediante la aplicación de tension y modificar también las propiedades elec- 
trónicas. Entonces, con el fin de describir los portadores de carga moviéndose en grafeno 
curvo, tomando en consideración además a los fonones flexurales; se ha obtenido una 
ecuación de Dirac generalizada. Para derivar dicha ecuación, se obtuvo primero una 
métrica apropiada y después se aplicó el principio de covarianza. La ecuación resultante 
contiene una corrección lineal, además de varios términos no lineales. Algunos de es- 
tos términos no lineales corresponden al desarrollo de Taylor del enfoque integral de 
traslape [14], mientras que otros, incluyendo la correccin lineal, son nuevos. Tales térmi- 
nos pueden dar lugar a efectos interesantes, como las resonancias entre los modos de 
fonones diferentes, los fonones debido a la flexión y el campo electromagnético externo. 
Esperamos que estos términos no lineales son importantes a temperaturas más altas. En 


contraste, tomando el caso particular, donde el grafeno es modelado por medio de una 
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forma tipo sinusoidal, se obtuvo un umbral para la propagación de electrones. Para los 
electrones por encima de este umbral, las ecuaciones correspondientes son sólo un con- 
junto de ecuaciones acopladas tipo Mathieu, lo que indica la posibilidad de resonancias 
paramétricas. Tales resonancias paramétricas pordrían dar lugar a grandes respuestas 
en el sistema. Esto es muy importante para establecer consideraciones de diseño de un 
dispositivo. No obstante, a partir de la relación de recurrencia, que surge de las ecua- 
ciones de Mathieu, pueden aparecer regiones de energías prohibidas [36]. El que exista 
regiones de energías prohibidas (GAP), podría permitirnos controlar y modificar el flujo 
de los electrones, y ser útil para construir compuertas electrónicas y diseñar transistores 
de grafeno. Esto es de gran importancia para las nuevas aplicaciones tecnológicas que 
se están llevando a cabo hoy en día. Cabe mencionar por último, que los resultados 


obtenidos se publicaron en la revista Physica B: Condensed Matter (ver Ref. [42]). 


Apéndice A 


Ecuación de Dirac 


Teniendo en mente el objetivo de resolver la ecuación de Dirac en un espacio plano, 
primero se realizará un análisis de dicha ecuación en 341 dimensiones. Posteriormente 
estos resultados serán trasladados al caso de 2+1 dimensiones y se comparará con la 


ecuación de amarre fuerte para el grafeno en los puntos de Dirac. 


A.l. Ecuación de Dirac en 3+1 dimensiones 


En 1927 Paul A. Dirac intentó postular una nueva ecuación para el movimiento de 
los electrones, que generalizara la ecuación de Schródinger, de modo tal que resultara 
compatible con la Relatividad Especial de Einstein (la Schródinger no era). Así que, para 
conciliar de manera consistente los postulados de la Relatividad con la interpretación 
probabilística de la teoría de Schródinger de la Mecánica Cuántica, Dirac optó por 


obtener de la ecuación de Schródinger 
ih— =HAV, (A.1) 


una ecuación invariante de Lorentz, que al ser lineal en la derivada temporal también 


debía ser lineal en las derivadas espaciales, y escribió la ecuación: 


HV =cd-pV + BmcV = O cop") Y + Bmc Y, (A.2) 


i=1 


donde c es la velocidad de la luz en el vacío, m la masa de la partícula, p el momento y 


a y B coeficientes a determinar. 
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En esta ecuación al imponer que se verifique la relación relativista 


PE=prime, (A.3) 
se tiene que 
8 | ] ze ] 
H=¿ y ¿(ar + ay,.04)pip* +me ) (048 + Ba)p' + Bm? (A.4) 
ik i=1 


donde ayajpip* = 5(oaj0jp ip? + ajop"p*) = B(ojayp ip? + a0pp*), pues i y k son 


índices mudos y p' y p" conmutan. Por tanto se deben verificar las condiciones, 


QiQAk + ARQ = 204 (A.5) 
o0B+ Ba; =0 (A.6) 
e =8=1 (A.7) 


Para cumplir estas condiciones los coeficientes q, y 6 deben ser matrices. Como condición 
adicional, las matrices deben ser hermitianas, para que el Hamiltoniano también lo sea. 


Ahora bien, usando el hecho de que a, = 8? = 1, obtenemos que los eigenvalores de 





estas matrices son +1. Además, tenemos que la traza de estas matrices es cero, esto se 


debe a la Ec. A.6, donde se encuentra que 
as = —BosP (A.8) 
Por las propiedades cíclicas de la traza se puede escribir 
Tra; = —TrBasB = —Trf*a; =0. (A.9) 


Ya que la traza es la suma de los eigenvalores, el número de eigenvalores positivos y 
negativos debe de ser igual, y como consecuencia a y Pf tienen que ser matrices con 
dimensión par. La mínima dimensión en que se cumple la relación 0/0 + 0x0; = 205% 


es 2 y la cumplen las matrices de Pauli 


ali) rm) en 


Sin embargo no existe ninguna matriz que cumpla 8? =1 y aB+Ba; = 0. Para cumplir 


estas dos condiciones es necesario ir a dimensión 4 y una representación particular de 


a ( 0), a=(; j . (A.11) 
(07 0 oO —1 


estas matrices es 
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donde las od; son las matrices de Pauli y las entradas con el símbolo 1 en la matriz /P 


representa matrices unitarias de 2x2. 


Así la nueva ecuación tiene 4 componentes y resulta 


Ov 


(CA -p+Bmc?) Y = No 


(A.12) 


donde el Hamiltoniano es 
H=cd-p+fme?. (A.13) 


El operador momento y la función de onda vienen dados por 


h > Wo 


(ñ,t) 

py | 24 (4.14) 
(ñ,t) 

(5,4) 


A.1.1. Forma covariante de la Ecuación de Dirac en 3+1 dimensiones 


Cabe señalar que la Ec. A.12 se puede escribir en forma covariante, es decir, si multipli- 


camos por £P a la izquierda de esta ecuación obtenemos: 
(cBa -p + me?) Y = mao (A.15) 
que puede reescribirse de la siguiente manera 
(cyp, — mec?) v =0, (A.16) 
(ihcy"9,, — mc?) Yv =0, (A.17) 
donde se han introducido las matrices de Dirac 
y = (8, Ba). (A.18) 


Las Ecs. A.16 y A.17 son llamadas comunmente como la ecuación de Dirac en su forma 


covariante. 


A continuación vamos a resolver la ecuación de Dirac para un partícula que se mueve 


en un espacio plano. 
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A.1.2. Partícula en movimiento 


Para describir a una partícula en movimiento, se propone una solución del tipo: 
Va) = HE TFDp(p). (A.19) 
Introduciendo esta solución en la Ec. A.16, obtenemos 


(cy"p, — mc?)p(p) =0. (A.20) 


Esta es la ecuación de Dirac en el espacio de momentos y 


se llama espinor de Dirac. 
Sustituyendo y“p, = Who — 7 -p en la Ec. A.20, llegamos a la siguiente expresión: 
E-m? -—có-p 
P op) Y _ o. (21) 
ap. —E-=mé x(p) 


5-5 cie ¡id (A.22) 


Px + 1Py —Pz 


con 


Este sistema posee 4 soluciones linealmente independientes: 


vale) = er = ur(pjeP?, (A.23) 





Ya(1) = rooms, (A.24) 
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3(1) = er? =us(pje Ps, (A.25) 


daa) =| Ent [| er? =us(pje >, (A.26) 


Interpretamos estas soluciones de tal manera que u¡ y uz corresponden al fermión (por 
ejemplo el electrón), mientras que uz y us corresponden al anti-fermión (por ejemplo el 
positrón). Las dos soluciones para cada caso representan a las dos direcciones de espín 


que puede tener la partícula. 


A.2. Ecuación de Dirac en 2+1 dimensiones 


Para obtener la ecuación de Dirac en 2+1 dimensiones, se debe eliminar la tercera 
componente espacial de la Ec. A.22 en 341 dimensiones; por lo que, sólo basta tomar 
ciertas matrices de 2x2 para satisfacer las condiciones A.5, A.6 y A.7. Estas matrices se 
pueden construir a partir de las matrices de Pauli, lo que da como resultado a que existan 


dos representaciones de las matrices y". Podemos escogerlas de la siguiente manera: 
P=0, Y =101, y =i02 (A.27) 
Y=o, y =i01, y =-i03. (A.28) 


Partiendo de la representación A.28 y las Ecs. A.21 y A.22, y proponiendo una solución 


de la forma A.19, obtenemos en el espacio momentos 


E-—m? 0 e 0 ¡Pa — Py ó pe (4.29) 
0 —E — m? ¡Py + Py 0 Xx 


Este sistema posee dos soluciones linealmente independientes: 








11 (x) = an ) e tp a ur(p)eP?, (A.30) 


E+mc? 
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—Py HF UPx 
Ya(x) = en ) e Pz = ua(pje*P, (A.31) 


La primera solución u¡ corresponde al fermión y la segunda solución us al anti-fermión. 


A.3. Ecuación de Dirac en grafeno plano 


Se sabe que las propiedades de transporte del grafeno, se centran entorno a los puntos 
K y K', llamados puntos de Dirac [14]. A partir de un estudio alrededor de ellos, se 
puede obtener la relación de dispersión, es decir, las bandas de energía; el cual son 
indispensables para saber si se trata de un metal o un semiconductor. Es preciso enfatizar 
que estos resultados surgen de resolver la ecuaciónn de Schródinger con la aproximación 
de amarre fuerte (tight binding) [43]. 


El modelo de amarre fuerte para el grafeno con interacciones a primeros vecinos, da una 


ecuación de Schródinger de la forma [43]: 


Hv(K)=E(K)V(K) y v(K)= eu ) l (A.32) 
x(K) 


donde el hamiltoniano es 


0 —at E ; 
H= E con a(K)=1+e Ka 4 ¿Kaz, (A.33) 
=axt 0 


La ecuación anterior pede reescribirse de la forma: 


E(K —at K 
(K) PK) =0. (A.34) 

=axt E(K) x(K) 
Lo anterior se vuelve un problema de eigenvalores, asimismo la relación de dispersión es 
obtenida en la forma estándar, haciendo el determinante cero. Por lo tanto, se tiene que 


la relación de dispersión es: 
E(K) = +ttla(K)| (A.35) 





Calculando la magnitud de a y usando las componentes (x, y) para K, obtenemos [43]: 








3 
DK By)= emp + 400 E, cos(EK,) + 4cos*(5 Ko), (A.36) 


donde a = yv3ag es la constante de la red y ay = 1,424 la distancia carbono-carbono en 


el grafeno. 
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Resulta interesante hacer un análisis más detallado de lo que sucede alrededor de los 
puntos de Dirac. Esto se puede hacer a partir de un desarrollo de Taylor, por ejemplo, 
del punto K. La función que vamos a desarrollar en serie en nuestro caso es a(K), el 


cual puede reescribirse como: 


ER 
a(K,, K,)=1+ 2e Ke cos(K). (A.37) 


En nuestro caso se toma las coordenadas del punto: 


27 27 


= a a 


Tenemos que la expansión en serie requerida alrededor de K + k es: 


.. (A.38) 


a(K + k) [iv 3aeTcos(Z), ae" sen(5) “(Rss lg) 


v3a 
2 





(úl — ky) (A.39) 


Si se toman nuevos ejes, de tal manera que el eje y se transforme en el eje x” y el eje x 


sea el eje y”, el hamiltoniano de la Ec. A.33 finalmente toma la forma: 


0 Billar — ik 
ya PE eN E (A.40) 
2 h(kg! + iky) 0 





donde up = v3at/2 = 1x10%m/s es la velocidad de Fermi. 


Notemos que el hamiltoniano de Dirac está dado por 
Hp =cd-p+ fBme?. (A.41) 


Si tomamos m = 0 y la representación de la ecuación de Dirac en 2+1 dimensiones, se 
tiene que: 


0 r=ib 
| ( añ (A.42) 
(Da! + iPy) 0 


Se observa qie la Ec. A.42 y la Ec. A.40 son equivalentes salvo que, para el grafeno la 
velocidad de los portadores de carga es la velocidad de Fermi up y que es 300 veces menor 
que la velocidad de la luz c. Por lo tanto, se puede estudiar al grafeno con la ecuación 
de Dirac ultrarelativista (como los fotones). Para el caso del grafeno los portadores de 


carga son llamados fermiones de Dirac sin masa. 


Apéndice B 


Ecuación covariante de Dirac en 


espacios curvos 


En esta sección comenzamos la búsqueda de las métricas resultantes y sus conexiones 
de Christoffel en grafeno corrugado. La superficie materializada por la hoja de grafeno 


puede ser descrito, por 


23= FE) (B.1) 


donde z, y son las coordenadas en la hoja de grafeno, y z es el plano de desplazamiento. 

Nótese que aquí sólo tenemos en cuenta la rama de fonones flexurales, ya que los fonones 

planos pueden ser descrito por la inclusión de un potencial vectorial [14], cuya solución 

es conocida [35, 36]. 

Como ahora la diferencial dz se convierte en una combinacin lineal de dx y dy, tenemos 
dz = dx + —dy (B.2) 

Por otro lado la métrica inducida en la hoja viene dada por la siguiente fórmula: 


ds? > gijdada?, (6,7 = 2,Y, z), (B.3) 


donde gi; es el tensor métrico. O bien, en una forma más explícita, 


ds? = da? + dy? + (Gea + qe ») ; (B.4) 


lo cual es una métrica en un espacio curvo de dos dimensiones, parametrizado por dos 


variables espaciales (x,y). Utilizando la última expresión se obtiene la forma explícita 
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de la métrica inducida, que es 


Ofy2  Ofof 
0 [ap da dy ) (B.5) 
of Of Ofy2 : 
so lA 


Esto puede ser escrito como 
Jik = 9jk + Ask con de k = L, Y, (B.6) 


donde gi; = diag(1,1) es la métrica plana en un espacio bidimensional, y h;, es la 
perturbación (se supone pequeña en comparación con 1) provocado por la ondulación 
de la hoja. La matriz B.5 es simétrica y real, por lo tanto puede ser diagonalizada por 
una transformación lineal apropiada. Así mismo la métrica inversa (contravariante) se 
puede encontrar fácilmente. Notemos que el determinante de la matriz correspondiente 


al tensor métrico covariante está como 


ña (y 50 2 2 
A (B.7) 
Ox Oy dy 


Por lo tanto, la matriz inversa, correspondiente a la métrica contravariante g/* es de la 








forma: 
p (9fy2 Of Of 
uE 14 (35) Ox Oy ) (B 8) 
of of p (9fy2 E 
Q dx dy 14 (ay) 


donde hemos utilizado la notación abreviada para el determinante, Q = (1 + (0, f)? + 
(9,f)?). La deformación hy se compone enteramente de términos cuadráticos, el cual 


va como el producto de las derivadas parciales espaciales de la función deformación f, 
hi0: 7059 (B.9) 
y es fácil comprobar que los símbolos de Christoffel correspondientes se reducen a 
Dir = 7 Om + Ordjm — Omán) = POS 05. (B.10) 


Ambas cantidades desaparecen cuando f = 0, y se convierte de nuevo en la métrica 
plana. La ecuación de Dirac en dos dimensiones debe ser ahora generalizada con el fin 
de incorporar el hecho de que la métrica en la superficie de restricción no es más plana, 


sino curvo. Ahora vamos a obtener la generalización covariante de la ecuación de Dirac. 


La ecuación que queremos producir puede ser escrita como la nueva deformación del 
hamiltoniano mecánico cuántico, actuando en un espinor Y de dos componentes de la 


siguiente manera: 


HU nr [UV + IV] U (B.11) 
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Aquí no sólo la métrica contravariante se deforma, sino también las matrices y deben 
ser modificados a fin de satisfacer las nuevas relaciones anticonmutativas con la métrica 
inducida, en lugar de la métrica plana; finalmente, Vj no sólo contiene la interacción 
electromagnética y el plano de fonones visualizado por el potencial vectorial, incluido de 
manera usual en la invarianza de norma, sino también los símbolos de Christoffel de la 
métrica Ji: 

V¡Y = (0, —eAj)U + Y. v, (B.12) 


donde los símbolos de Christoffel se definen como de costumbre, por medio de la métrica 


modificada: 
e Ls Ñ - ) 
Da E UN [OjImk + Or Ijm — OmIjk] (B.13) 


y m0 es el valor matricial del tensor antisimétrico definido por medio de las matrices- 


gamas modificadas: 


se o e 


Este término a menudo se llama ”conexión espinorial” [37, 38]. 


Estamos buscando a dos generadores deformado” del álgebra de Clifford Y, y Yy, que 


satisfagan 


Of 0 
Ya y + VW Ye = 2 ol. (B.15) 


Con el fin de hacer esto, hay que introducir la tercera matriz de Pauli, ya que la defor- 
mación de la lámina se sale fuera del plano bidimensional (x, y). Por otro lado, el álgebra 
de Clifford sin deformar que satisface las relaciones anticonmutativas en una red plana 


en el espacio de tres dimensiones, se define como sigue: 
ViVe + Ve Y =29jr1 con gjx=diag(1,1,1), ¿,k=1,2,3, (B.16) 
que pueden ser generados por tres matrices de Pauli como sigue: 
VU Y =Uy, Y3=0z (B.17) 


con 


y además 
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El ansatz para las dos deformaciones en el espacio donde viven las matrices-y es simple 
si ponemos 


Ye = 0 H000 Ty = 0 FO (B.20) 
Entonces los coeficientes de a y de b deben ser, como es fácil de comprobar, 


01 


de modo que 


0 Of 


Ahora tenemos que producir sus contrapartes contravariantes que aparecen en la ecuación 


Ya = Og + Oz, VYy = 0y + 


de Dirac (ver ápendice A). Tenemos 


Y = GIA A GU, PA GA A GUA 45048) 
que da explícitamente 

5 = ¿[1+(0, 0, +(0,4)0, —(0:0,P)04). 

ai 5 [1 + (9,170, +(0,$)0 — (0,$)(0,$)0%] - (B.24) 


Recordando que Q = (1+ (09, f)? +(0,.f)?), podemos sumar y restar en los numeradores 
de la fórmula anterior, respectivamente, los siguientes términos: (9, f)?0, en la primera, 
y (0 PDA en la segunda; lo que nos permitirá separar las matrices no deformadas 0. y 
0y, y los términos originales de la deformación que contienen derivadas espaciales de f. 


Esto da el siguiente resultado: 

ES pa GS e => 1 an ES 2 

Hr pWO y ot Vy + (9rad] -V)— y gradf -g)gradf -V), (B.25) 
donde los vectores y sus productos escalares son de dos dimensiones, es decir, 


gradf =(0+f,0%f), 5=[0%,0y), V= (Va, Vy] 


de modo que 
gradf -?=0,f0% + Ou iO Ele (B.26) 


Vemos que ya en los numeradores tenemos no sólo términos lineales, sino también térmi- 
nos cuadráticos, a pesar de la presencia de términos cuadráticos en el denominador 
(contenida en el factor de normalizacin de (2). Si decidimos mantener sólo términos lin- 


eales, entonces el hamiltoniano modificado debería contener sólo un término adicional 
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proporcional a la matriz 0: 
Bin - Mo +85 =0*"V¿ +09V, +0: (gradf - V). (B.27) 


Tenga en cuenta que también el operador diferencial se toma en su forma primaria, 
debido a que los coeficientes de la conexión sólo contienen expresiones cuadráticas en 
las derivadas de f. Observe que esta ecuación es similar a la obtenida en el desarrollo de 
Taylor en el enfoque de amarre fuerte. Sin embargo, si optamos por mantener a todos los 
términos de segundo grado, entonces debemos también tener en cuenta los coeficientes 
de Christoffel en V. Es fácil comprobar la siguiente forma explícita de los símbolos 
de Christoffel; conservando sólo las expresiones de segundo orden significa que podemos 


utilizar la fórmula simplificada en la que 37 se sustituye por gY = 9. Tenemos entonces 
E-= 00d Cy = Pia = 0:f0%,F, Oy = O f0%,f, 


O E OO de 10 O (B.28) 


En derivadas covariantes, estos coeficientes se contraen con g/* y las matrices antisimtri- 


cas Y tm» 


Estamos utilizando las matrices no deformadas >, jk Y nO los deformados $ jk» Porque 
vamos a mantener las expresiones hasta la segunda potencia de las derivadas de f, que 
ya están contenidos en los símbolos de Christoffel. Esto, teniendo en cuenta que sólo 


los términos diagonales en g?* desaparecen y son iguales a uno, conduce al resultado 


siguiente. Para j = 1 tenemos 


Dig Dan +39 Y ay + Dtag7 Y ya + Digg Y py 


y para j = y se obtiene una expresión similar: 


A, x Y TL y 
Dag Y a + Py 9 Y ay + ya Y ya + Di 9 Y y 
Los componentes del tensor métrico no se desvanecen y son iguales a uno, mientras que 


las matrices » son antisimétricas en sus dos índices más bajos, así que lo que queda es 


sólo 


Ey y TY - Das y yr — (0-f02,f > Oy F Oz f) Y TY 
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Te, Y ay HT Y Y yo = (04108, $ — 00482,$) Y y (B.29) 


Estos son los únicos términos restantes a ser incluidos en las derivadas covariantes. Estos 


se escriben como sigue: 


Va = (07 — € Ag) + (0 0%yf — 014 0%2f) Y ay) 
Vy = (9, — e Ay) + (0,40%, f — 040%, 4) Y yo: (B.30) 


En el hamiltoniano, ellos parecen multiplicarse por la izquierda por las correspondientes 
matrices y, pero aquí, la evaluación de los términos procedentes de los coeficientes de la 
conexión de Christoffel, ya de segundo grado en la deformación f, se puede conservar 
sólo en su versión deformada, que en nuestro caso son las dos matrices de Pauli o” y 0”, 
por lo que la parte del hamiltoniano deformado manteniendo la conexión de Christoffel 


es 


NV UN av, = o* (0, E eAz) 4 jos => Pron) EN zY 
+0Y (0, — Ay) + 09 (0,4 0%0f — 059%) Y ya (B.31) 


Teniendo en cuenta el hecho de que 


1 
lu) =1, (a) =1l, 0501 =-9y00, y 0% P. = 70 (B.32) 


2 . ., 2 
sumando y restando 0, f0;,, f con la primera expresión y 0, f0;,, f con la segunda, obten- 
emos la siguiente forma invariante de los términos adicionales inducidos por la conexión 


de Christoffel: 


1 
gor (Of 0% t +0,8%y $ — 0450 — 040%) 


+ (Op$0%,$ + 005023 — O0f0%af — 0:1%,S) (B.33) 


Esto a su vez puede ser escrita en una ecuación más compacta (y elegante) de la manera 
siguiente: 


37 gradlo.1 Y + (0,1710 -gradfjAs (B34) 


1 
4 


con 
Af=0O%0f + Oiyf 


Vale la pena señalar que la Ec. B.34 desaparece cuando f tiene una forma geométri- 
ca sinusoidal, pero es diferente de cero tan pronto como no exista una superposición 
de estas expresiones, por ejemplo, para una onda estacionaria. Ahora somos capaces 


de escribir el hamiltoniano total de un electrón en una hoja, teniendo en cuenta que 
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movimientos propios de la hoja descrita por la deformación del plano horizontal dada 


por z = f(x,y,t), hasta segundo orden (términos de segundo grado en las derivadas de 


pr 
Aa goalarads) . ¿(oradf rar) 


+7 -gradllo,1) +(0,171- 10 -gradf]A7 (B.35) 


1 
4 
El factor de normalización 1/Q en frente de las dos primeras contribuciones se puede 
ajustar a 1, ya que contiene los cuadrados de las derivadas de f, y si se desarrolla, va a 
crear condiciones de orden 3 y 4 cuando se multiplica por los términos que se encuentran 
a su derecha. El hamiltoniano de la mecánica cuántica se obtiene de esta expresión 


multiplicándolo por —¿ y vp. Vamos a definir el operador de momento generalizado, 


> 


ll =P-eA4= —ihgrad —- eA, 


donde Á es un potencial vectorial que describe un campo electromagnético [18] y [19] o 
planos de fonones longitudinales y transversales [17] . Entonces, si hacemos grad — V, 


el hamiltoniano toma la forma: 


=p Y, ) e (Os f)%z + (04 f)%y 0 ] 
Ú (Oz f)Tz = (Oy f)Ty 


0 p(ÚS GNU 4) — ihr [ho VOLS) +(9,87] - 415 - VIAS. (B.36) 


Esta última ecuación es la ecuación de Dirac generalizada. Es interesante observar que el 
hamiltoniano que aquí se deriva, la corrección a primer orden, los términos son similares 
a los obtenidos a partir de un enfoque de potencial vectorial [14]. Sin embargo, la Ec. 


B.36 contiene más términos. 
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